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Abstract 

Large-scale motions, important in turbulent shear flows, are frequently attributed to the inter- 
action of structures at smaller scale. Here we show that, in a turbulent channel at Re T ~ 550, 
large-scale motions can self-sustain even when smaller-scale structures populating the near-wall 
and logarithmic regions are artificially quenched. This large-scale self-sustained mechanism is not 
active in periodic boxes of width smaller than L z ~ 1.5/j or length shorter than L x ~ 3h which 
correspond well to the most energetic large scales observed in the turbulent channel. 



Understanding the mechanisms by which turbulence self-sustains in wall bounded tur- 
bulent shear flows is still an open challenge that is relevant for applications ranging from 
atmospheric dynamics to the design of low drag vehicles. One of the most robust features of 
these flows is the presence of streamwise streaks, i.e. narrow streamwise regions where the 
streamwise velocity is larger or smaller than the mean velocity at the same distance from the 
wall. The existence of streaks in the buffer layer with a characteristic average spanwise scale 

n n 

A+ ~ 100, in wall units, has been revealed by early experimental observations [l|, [2]. Nu- 
merical experiments furthermore reveal that turbulence is not sustained in periodic domains 

n 

smaller than this spacing |3J. There is now large consensus that these near- wall structures 
are generated via a self-sustained process 14, 5 1 based on the streaks amplification from the 

n on 

vortices via the lift-up effect [6J followed by the breakdown of the streaks p, |8[ and the 
regeneration of the vortices. This buffer-layer process survives when turbulent structures 
are artificially removed from the outer layers Jq]. 

Coherent structures, however, exist at larger scales. It has long been known that stream- 
wise velocity correlations are important up to lengths of the order of the outer length scale 
h of the flow (e.g. the channel half-width, the pipe radius or the boundary layer thickness), 



and have been related to the presence of 'bulges' or 'large-scale motions' 



10j. Recently, it 



has been also realized that coherent streaks exist at even larger scale with typical spanwise 



and streamwise scales \ z ~ 0(h) and X x ~ O(10h) respectively ll|-15|. These very large- 
scale streaks, also referred to as very large-scale motions or global modes or superstructures, 



are important because they carry a very significant 



and of the turbulent Reynolds stress 



12 



13 



15 



Taction of the turbulent kinetic energy 



Ipl contradicting the early view that the 



17|. 



18|. 



motions at very large scale are essentially inactive 

The fundamental question in which we are interested here is: What is the origin of the 
streaky motions at very large scale? It has been conjectured that these motions could result 
from the concatenation of large-scale motions 13j, ll6[. These large-scale motions are in t 
turn made of the aggregation of a huge number of hairpin vortices of smaller scale 12L 
Even if from a different perspective, further numerical experiments 19| support the idea 
that very large scale motions can be directly forced and interact with near-wall structures, 
at much smaller scale, via a co-supporting cycle. In this mechanism, while the near-wall 
cycle would continue to exist even in the absence of large scale structures, the reverse would 
not happen. 
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The current wisdom is therefore that large- and very large-scale motions would not ex- 
ist in the absence of smaller scale structures. An alternative way of thought is however 
emerging from recent findings that extend recent hydrodynamic stability methodologies to 
fully developed turbulent flows. It has recently been theoretically predicted that very large 
scale streamwise streaks can be amplified by a coherent lift-up effect which is able to extract 
energy from the mean flow at very large scale, without the mediation of coherent structures 



at smaller scale 



20M24J|. The existence of the coherent lift-up has been confirmed in ex- 



periments where very large scale coherent streaks have been artificially induced by forcing 



coherent streamwise vortices 



25 



5l-l27l|. It is therefore tempting to conjecture that the very 
large scale motions could result from a self-sustained mechanism bearing some similarity to 
the processes observed in the buffer-layer and in transitional flows (see e.g. the discussions 



m 



13 



22 



281]) . The goal of the present study is to verify if this conjecture is true. 



We proceed by designing a 'conceptual' numerical experiment where we gradually remove 
the smaller scales from the flow while still taking into account the associated dissipation. The 
ideal tool to perform such a kind of experiment is the large eddy simulation (LES), which has 
also the advantage of keeping manageable the computations in very large domains at high 
Reynolds numbers. In LES the very small scale motions are not resolved by the grid and 
the residual stress tensor associated with these small scale filtered motions is modelled with 
an eddy viscosity ry — r kk 5ij/3 = —2v T Sij, where is the symmetric part of the filtered 
velocity gradient. From the adopted perspective it is important to avoid any backward 
spectral energy transfer from the smaller to the larger scales. Therefore, we consider the 
'static' Smagorinsky model (Smagorinsky, 1963) where the eddy viscosity is simply modeled 
as v T = D(Cs&) 2 S where, S = (2SijSij) l l 2 is a measure of the local shear and C5A is the 
local mixing length, given by the product of the filtering scale A to the Smagorinsky constant 
Cs- The van Driest damping function D = 1 — e~( y+ / A+ ^ 3 is also used to enforce a physical 



behavior near the wall 



or the further details on properties of the Smagorinsky model, the 



reader may refer e.g. to 29 



30J]). Dynamic Smagorinsky models, that are known to perform 



better for engineering applications, are not used here because we want to avoid backscatter. 
The simulations are performed using the public domain code diablo [31J which implements 
the fractional-step method based on a semi-implicit time integration scheme and a mixed 
finite-difference and Fourier discretization in space. The computational domain is chosen 
large enough (L z = 6h) to contain three to four large scale streaks in the spanwise direction 



3 




FIG. 1: Comparison of mean statistics of the reference case (C5 = 0.05) with data from the direct 



numerical simulations in Ref. 



331 ] at the same Re T = 550: (a) mean streamwise velocity profile 



U(y + ); (b) turbulence intensity profiles u rms (y + ), v rms (y + ) and w rms (y + ). 

and long enough (L x = 5Qh) to contain the longest observed streamwise wavelengths of the 
very large scale motions. The grid size of 384 x 65 x 96 points in respectively the streamwise 
(x), wall- normal (y) and spanwise (z) direction is set-up so as to resolve the near- wall cycle 



properly in the reference simulations 32]. A grid stretching is applied in the wall- normal 
direction with the first grid point at y + = 1.8 and a maximum spacing A+ = 40 at the 
channel center. 

The simulations are performed by enforcing a constant mass flux corresponding to the 
constant Re m = 20133, where Re m = U m 2h/u is the Reynolds number based on the bulk 
velocity U m = (l/2h) j\ U dy, the kinematic viscosity v and the channel half-height h. The 
reference case is first considered by choosing the Smagorinsky constant C$ = 0.05 which was 



shown to provide the best accordance with statistics of direct numerical simulations 29]. 
The reference solution exhibits an average friction-based Reynolds number Re T = 550, where 
Re T = hu T /v and u T is the friction velocity. The mean streamwise velocity profile and the 
turbulence intensities agree well with the those issued from direct numerical simulations at 



the same Re T [33j as reported in Fig. [TJ A snapshot of the streamwise velocity fluctuation in 
the reference case is reported in Fig. Eta). The structure of this field is very similar to the one 
obtained from direct numerical simulations, with large-scale streaks embedded in clouds of 



smaller scale structures (see e.g. 28]). The premultiplied spanwise and streamwise velocity 
spectra k z E uu {k z ) and k x E uu {k x ) corresponding to the reference case, reported in Fig. Eta) 
and Etc) respectively, are also in qualitatively good accordance with those issued from DNS 



at similar Re T Uj. In particular, in the premultiplied spectra, the scales of the buffer-layer 



cycle A+ ~ 100 and A+ ~ 600 are observed for y + < 30 , while in the outer region the peaks 
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FIG. 2: (Color online). Isosurface (u + = —2) of the instantaneous streamwise velocity fluctuation 
for: the reference simulation with Cs = 0.05 (a) and for simulations performed with increasing 
artificial dissipation of the small-scales with Cs = 0.1 (b), Cs = 0.2 (c) and Cs = 0.3 (d). 

corresponding to the large-scale motions are observed at X z /h m 1.5 ~ 2 and X x /h ~ 3 ~ 4 
with long tails in the range of X z /h > 3 ~ 4 that represent the very-large-scale motions 



16 



34 



35|. 



To test if the large-scale structures can survive even in the absence of smaller scales 
we progressively damp the small-scale structures by simply repeating the simulations with 
larger values of the Smagorinsky constant Cs which only affects the diffusion term and not 
the nonlinear advective terms. In these additional simulations the Reynolds number is kept 
at its reference value Re m = 20133. The friction-based Reynolds number is not strongly 
affected by the change in Cs- Re T = 494 is found for Cs = 0.1, Re T = 518 for Cs = 0.2 and 
Re T = 560 for C s = 0.3. 

For increasing Cs values, the small-scale coherent motions are increasingly damped, as 
can be seen from the instantaneous streamwise velocity fluctuation fields reported in Fig. [2] 
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Only large- and very large-scale streaky motions survive for Cs ~ 0.3 (Fig. [2$) and also these 
motions are quenched for Cs = 0.4. The effects observed on the instantaneous velocity fields 
are confirmed by the analysis of the premultiplied spectra reported in Fig. [3J For Cg = 0.1, 
the buffer-layer peak A 2 + ~ 100 in the spanwise premultiplied spectrum k z E uu (X z ) has been 
quenched, meaning that the near wall cycle has been suppressed. However, structures with 
almost constant X z /y survive in the log-layer for Cs = 0.1, as seen in Fig. [3^c) . For larger 
Cs, also these structures are quenched. Only one peak survives with X z w 1.5h for Cs = 0.2 
and X z = 2h for Cs = 0.3 with X z corresponding well to the scale of large- and very-large 
scale motions (A 2 ~ 1.5 — 2h) apparent in Fig. [Hfr and observed in numerical simulations 
and experiments 

HQ. 

The analysis of the streamwise premultiplied spectra k x E uu (X x ), confirms this scenario. 
For Cs = 0.1 the X x + ~ 600 peak has been replaced by much longer structures (X x ~ 7h). In 
the channel center, structures with X x ~ 3h dominate like in the reference simulation. The 
typical streamwise scales in the inner layer, no more hidden by the near-wall cycle, are now 
longer than in the outer layer. These scales grow longer when Cs is further increased (see 
Fig. |3](d),(f),(h)). The streamwise scales surviving for Cs = 0.3 are larger by a factor of two 
when compared to the 'natural' large-scale and very large-scale motions, which can probably 
be attributed to the very important artificial dissipation enforced at this Cs- The analysis of 
the wall-normal structure of the modes associated with these scales, reveals that the modes 
associated to longer X x , when compared to shorter X x ones: (a) have turbulence intensity 
peaks nearer to the wall, (b) penetrate more to the wall preserving a larger turbulence 
intensity, (c) have a ratio of the wall-normal to the streamwise turbulence intensity which is 
smaller. These structures remain essentially outer-layer structures even for the longest X x . 

These findings, especially for the case Cs = 0.1, are in accordance with recent numerical 
experiments where it was shown that lar ge s cale motions in the outer layers are almost 
independent of the buffer-layer processes 36j, |37|]. Even more importantly, these results 
prove that large- and very-large scale motions in the turbulent channel do not necessarily 
rely on the existence of smaller scale structures like buffer-layer streaks and vortices or 
hairpin vortices. This is a strong indication that an independent self-sustained process at 
large scale exists in the turbulent channel and probably in the other turbulent canonical shear 



flows such as the pipe and the boundary layer, as was previously conjectured 



13 



20 



22]. 



The self-sustained process at large scale is associated with coherent large scale streaks 
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FIG. 3: Spanwise premultiplied power spectrum k z E uu (X z ) [(a), (c), (e) and (g)] and streamwise 
premultiplied power spectrum k x E uu (X x ) [(b), (d), (f) and (h)] for respectively the reference simula- 
tion with Cs = 0.05 [(a) and (b)] and for the cases Cs = 0.1 [(c) and (d)], Cs = 0.2 [(e) and (f)] and 
Cs = 0.3 [(g) and (h)]. The premultiplied spectra are extracted in the inner- (y + = 16,30,70,108 
dashed lines) and in the outer layer (y/h = 0.38, 0.65, 1 solid lines). 

that undergo sinuous oscillations and breakdown in an apparently random way. This phe- 
nomenology is similar to what is observed in minimal box realizations of the buffer layer 
self-sustained process and in the subcritical laminar-turbulent transition in the channel 
0, 0, Q, 38 ] . In these situations, the process can self-sustain only in periodic boxes suffi- 
ciently large for the streaks to grow by lift-up effect and sufficiently long to allow the sinuous 
streaks oscillations to develop. We have therefore repeated the simulations for Cs = 0.3 in 
smaller horizontally-periodic boxes. A constant spatial resolution (Ax, Az) is used so as to 
maintain the LES filter width, unchanged. A detailed analysis of numerous box combina- 
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tions reveals that the large-scale process is sustained only if the streamwise and spanwise 
box sizes are larger than the minimal values L x>m i n « 3h and L Z)Tn i n fa 1.5h respectively. 
These values are in very good agreement with the spatial length-scales of most energetic 



outer layer motions 
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15 



20 



341 ] . This similarity of the minimal box sizes to the char- 



acteristic scales of the most energetic structures observed in the 'natural' turbulent flows is 
also observed for buffer-layer structures j^], The minimal spanwise size L z ^ m i n is also well 
in the range of scales where streaks are strongly amplified by the coherent lift-up effect in 
turbulent channels 



20 



211. 



The analysis of the power input and dissipation I(t) and D(t) in the minimal box L XtTnin x 
L z , m in for Cs = 0.3 reveals pseudo-periodic oscillations with period ~ 13 — 15/i/L^ associated 
with sinuous oscillations of the low speed streaks in the whole channel. These oscillations 
are typically followed, in an intermittent way, by a large excursion in the drag and the 
dissipation that is very similar to the 'bursting' observed in the buffer-layer process \3\. 

From the minimal box results we also conclude that the existence of motions at very large 
scale (X x > 3h) is not necessary for the survival of the motions at large scale (X x ~ 3h). 
Whether the motions with X x > 3h are active or are passively forced by the process at 
X x ~ 3h, is still an open question that is under current scrutiny. 
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part of this work was done. 



* Electronic address: yongyun@ladhyx.polytechnique.fr 



t Electronic address: Icarlo . cossu@imft irl 
[1] S. J. Kline, W. C. Reynolds, F. A. Schraub, and P. W. Runstadler, J. Fluid Mech. 30, 741 
(1967). 

[2] J. R. Smith and S. P. Metzler, J. Fluid Mech. 129, 27 (1983). 

[3] J. Jimenez and P. Moin, J. Fluid Mech. 225, 213 (1991). 

[4] J. Hamilton, J. Kim, and F. Waleffe, J. Fluid Mech 287, 317 (1995). 

[5] F. Waleffe, Phys. Rev. Lett. 81, 4140 (1998). 

[6] M. T. Landahl, J. Fluid Mech. 212, 593 (1990). 



S 



[7] F. Waleffe, Stud. Appl. Math. 95, 319 (1995). 

[8] W. Schoppa and F. Hussain, J. Fluid Mech. 453, 57 (2002). 

[9] J. Jimenez and A. Pinelli, J. Fluid Mech. 389, 335 (1999). 

[10] L. S. G. Kovasznay, V. Kibens, and R. F. Blackwelder, J. Fluid Mech. 41, 283 (1970). 

[11] J. Jimenez, Annual Research Briefs, Center for Turbulence Research, Stanford University 
(1998). 

[12] C. D. Tomkins and R. J. Adrian, J. Fluid Mech. 490, 37 (2003). 

[13] M. Guala, S. E. Hommema, and R. J. Adrian, J. Fluid Mech. 554, 521 (2006). 

[14] N. Hutchins and I. Marusic, J. Fluid Mech. 579, 1 (2007). 

[15] B. J. Balakumar and R. J. Adrian, Phil. Trans. R. Soc. A 365, 665 (2007). 

[16] K. C. Kim and R. Adrian, Phys. Fluids 11, 417 (1999). 

[17] A. Townsend, The structure of turbulent shear flow (Cambridge U. Press, 1976), second edi- 
tion. 

[18] X. Wu and P. Moin, J. Fluid Mech. 630, 5 (2009). 

[19] S. Toh and T. Itano, J. Fluid Mech. 524, 249 (2005). 

[20] J. C. del Alamo and J. Jimenez, J. Fluid Mech. 559, 205 (2006). 

[21] G. Pujals, M. Garcia- Villalba, C. Cossu, and S. Depardon, Phys. Fluids 21, 015109 (2009). 

[22] C. Cossu, G. Pujals, and S. Depardon, J. Fluid Mech. 619, 79 (2009). 

[23] Y. Hwang and C. Cossu, J. Fluid Mech. 643, 333 (2010). 

[24] Y. Hwang and C. Cossu, Submitted to J. Fluid Mech. (2010). 

[25] O. Kitoh and M. Umeki, Phys. Fluids 20, 025107 (2008). 

[26] G. Pujals, C. Cossu, and S. Depardon, in Sixth Syrup, on Turbulence and Shear Flow Phe- 
nomena (Seoul Nat. University, Seoul Korea, 2009). 

[27] G. Pujals, C. Cossu, and S. Depardon, Submitted to J. of Turbulence (2010). 

[28] J. Jimenez, Rev. R. Acad. Cien. Serie A Mat. 101, 187 (2007). 

[29] C. Hartel and L. Kleiser, J. Fluid Mech 356, 327 (1998). 

[30] C. Meneveau and J. Katz, Annu. Rev. Fluid. Mech. 32, 1 (2000). 

[31] T. R. Bewley, P. Moin, and R. Temam, J. Fluid Mech. 447, 179 (2001). 

[32] T. A. Zang, Phil. Trans. R. Soc. Lond. A 336, 95 (1991). 

[33] J. del Alamo and J. Jimenez, Phys. Fluids 15, L41 (2003). 

[34] J. C. del Alamo, J. Jimenez, P. Zandonade, and R. D. Moser, J. Fluid Mech. 500, 135 (2004). 

9 



[35] J. P. Monty, N. Hutchins, H. C. H. NG, I. Marusic, and M. S. Chong, J. Fluid Mech. 632, 
431442 (2009). 

[36] O. Flores and J. Jimenez, J. Fluid Mech. 566, 357 (2006). 
[37] O. Flores, J. Jimenez, and J. del Alamo, J. Fluid Mech. 591, 145 (2007). 
[38] S. C. Reddy, P. J. Schmid, J. S. Baggett, and D. S. Henningson, J. Fluid Mech. 365, 269 
(1998). 



10 



